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The quasiclassically exact solution for the second inertial parameter B is 
found in a self-consistent way. It is shown that superdeformation and nonuni- 
form pairing arising from the rotation induced pair density significantly reduce 
this inertial parameter. The new signature for the transition from pairing to 
normal phase is suggested in terms of the variation B/A versus spin. Exper- 
imental data indicate the existence of such transition in the three SD mass 
regions. 

PACS numbers: 21.10.Re, 21.60.Fw, 23.20.Lv, 27.60.+j 

One of the amazing features of superdeformed (SD) rotational bands is the ex- 
treme regularity of their rotational spectra: a SD nucleus is the best quantum 
rotor known in nature. In spite of the fact that numerous theoretical calculations 
successfully reproduce the measured intraband 7-ray energies (see e.g. [1-4]) the 
underlying microscopic mechanism of this phenomenon is still not well understood. 
To explain the SD rotational spectrum regularity we use in the present paper the 
parameterization of its energy by the three term formula 
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E(I) = E + AI(I+l) + BI 2 (I + l) 2 , (1) 
which is valid for an axially symmetric deformed nucleus with K — 0. The inertial 



parameters A = h /2Q 



is the kinematic moment of inertia) and B are 



the objects of our investigation. They are determined by the transition energies 
E 7 (I) = E(I + 2) - E(I) as follows 
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The coefficient B(I) characterizes the nonadiabatic properties of a band and is very 
sensitive to its internal structure. In particular, it realizes the relationship of kine- 
matic and dynamic (Q^ 2 )) moments of inertia. Using the well known expressions for 
these values || and the formulas (|2|) we get 
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The ratio B/A determines the convergence radius || of the two parameter formula 
(HD, which is of order 100 for the bands in the 80 and 150 mass regions and 40 in 
the 130 and 190 ones. Faster convergence is obtained with Harris formula 

E(u) = E + X -au? + (4) 
which is based on the fourth-order cranking expansion 

« = -5p(4p {1) ), (3 = -^Sp(l xP ^), (5) 

where the nth correction to the nucleus density matrix, £ x is the single- 

particle angular momentum projection on the rotational axis x perpendicular to the 
symmetry axis z, and u is the rotational frequency. It follows from Eqs. ([!]) and (S) 
that 

a = h 2 /2A, (3 = -h A B/AA 4 . (6) 

For simplicity we will deal with the parameter j3. 

The problem of the microscopic calculation of the parameter B for normal de- 
formed (ND) nuclei has attracted considerable attention. Its value is formed mainly 
by four effects: vibration-rotation interaction, centrifugal stretching, perturbation 
of the quasiparticle motion, and attenuation of pair correlation by the Coriolis force 
(Mottelson-Valatin effect). As it has already been shown in the first attempts at 
obtaining B |J7|, the latter two effects are dominant for well deformed nuclei. 
Unfortunately the results of these and many other works cannot be used for the 
superdeformation. The formulas of Ref. || have been obtained in the limit of the 
monopole pairing interaction (the uniform pairing), which is not adequate at SD 
shapes as shown in Ref. . In the more sophisticated work |7| the gauge invariant 
pairing interaction allows to study the effect of nonuniform pairing. However this 
approach is also limited because it neglects the coupling between major shells. Thus, 
despite a number of publications on the subject the correct cranking selfconsistent 
solution for the B coefficient has not been found. 

We used the quasiclassical method of Ref. to derive the following expression 
for the (3 parameter in the superfluid phase 

& = - 4^2 H £ 12^23^34^41 F 12 , ^23 , ^34 , X 41 ) 6 (f- 1~ £ F ) , (7) 

where the summation indices i=l, 2,3,4 refer to the single-particle states i of the 
nonrotating mean field with the energy £j. The 5-function means that the summation 
over the states 1 is restricted by a small interval at the Fermi energy Bp [0. The 
dimensionless values xu> = (e* — £j/)/2A, where A is the state independent pairing 
gap at w = 0, correspond to energy differences between states permitted by the 
selection rules for the matrix element of i x . The function F depending on these 
values may be written as follows 

3 1 

F=^2P k G 12 +J2 p k H i3-8Dlx 12 X23X U x 41 h(xi 3 ), (8) 

k=0 k=0 
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where the permutation operators in the space of four indices i, i modA = i, 

Pk%i,i' = x i+k,i'+ki (9) 

are used to simlify the formulas. The expressions for Gxi and H13 involve the well 
known functions ]10| 



h(x) = (l + x 2 )g(x), g(x)= "T^ , (10) 



and have the form: 



G 12 = 9 ^ l2> f(l-D l xl 2 )\-l-xl 2 -X23X4l+D 1 [xl 3 (l-X 12 X24)+xl 1 (l+X 1 2X l3 ) 



^23^41^13^24 

+ Dlx2 3 x u x 41 (x23+x 41 )Dfx 12 x 2 23 x u xl 1 + £>i(x 34 - D 1 X 12 X2zX il ){x M 

+ X 12 X 13 X 2 4 ~ -D1X12X23X41)} , 
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Eq. (0) multiplied by a; 3 is the third order cranking correction to the total an- 
gular momentum of the neutron or proton system. Its derivation will be described 
in a forthcoming paper. We want now to emphasize that Eq. (|7]) represents the 
first theoretically correct expression for the high order effect of the Coriolis-pairing 
interaction at fixed deformation. The result is obtained by taking into account the 
effect of rotation on the Cooper pairs in the gauge invariant form. This effect is 
described by the first and the second corrections to the pairing energy 

A«W = -^i4 A«(r) = ?^4, (12) 

where Di and D 2 are the amplitudes of the nonuniform pairing fields, which are 
found in a self-consistent way. Let us note that is the function of the space 
coordinates because M x = ydU/dz — zdU/dy, where U is the potential of a mean 
field. It is seen that for the oscillator potential l x ~ yz (the Y 2 \ pairing). The 
coordinate dependent pairing field is crucial for conservation of a nucleon current. 
The theory incorporating the nonuniform pairing allows also to consider the different 
limiting cases for the inertial parameters, which make possible the study of an 
interplay between rotation, pairing correlations, and mean field deformation in a SD 
band. 

In order to consider this problem quantitatively we will use the axially deformed 
oscillator potential with the frequencies uj x = uo y and uo z on the corresponding axes. 
In this model the matrix element l\ 2 is non-zero for the two types of transitions: 
ii) transitions inside a single oscillator shell (close transitions), for which Xi2 = ±^i; 

transitions over a shell (distant transitions) with Xi2 = ±^2- The quantities v\ 
and v 2 are the well known parameters involved in the moment of inertia JT0]|: 



Ul ' 2 ~ 2 A -2^' (13) 
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where hu = AlA^ 1 ^ 3 MeV. Here and later we use the axis or frequency ratio k = 
c/a = oj x /uj z and the volume conservation condition. Both of the values v\ and 
z/ 2 are large for superdeformation. For the fixed state 1, there are 36 different 
combinations of these base transitions in the sum of Eq. (0). The summation over 
the states 1 is performed in the Thomas-Fermi approximation. The final expression 
for the parameter f} s in the oscillator potential is 

where R = l.2A 1 ' 3 fm is the radius of the sphere, which volume is equal to that of 
the spheroid with the half-axes a < c, M is the nucleon mass, and A is the number of 
nucleons. The function $ along with its limiting cases is shown in Fig. 1. It is seen 
that nonuniform pairing reduces substantially the parameter (3 S in agreement with 
the estimation of Ref . || . On the other hand, the contribution of distant transitions 
is minor for small £. Nevertheless the later are necessary to obtain the hydrodynamic 
limit (see below). Since $ ~ 1 for a reasonable pairing gap, A ~0.5MeV, the order 
of the value (3 S is h A (A/ Sp) 3 ■ This, along with the estimation A ~ epA~ 5 ^ 3 , gives 
B/A ~ A~ 2 , which overestimates the minimal value of this ratio in all the SD mass 
regions. Thus a small A and nonuniform pairing does not solve the problem of the 
SD band regularity. 

Let us consider first the limiting case A = 0. The right part of Eq. (0) vanishes for 
noncorrelated nucleons. This result is the artifact of the quasiclassical approximation 
used in Eq. ([?]). The correct expression for the f3 parameter in the normal phase, 
obtained with the limiting values of the Bogolubov amplitudes (ui = 0, Vi = 1 for 
rii = 1 and Ui = 1, Vi = for n, = 0, where rii is the nucleon occupation numbers), 
has the form 

k=0 £ 12 £ 13 £ 14 

The odd function of the differences eu> = e, — leads to the cancelation of the main 
terms in sum (|l^) that decreases substantially the value of /3 n , (3 n ~ h A A 7 ^/e i F . Note 
that the centrifugal stretching effect has the same order /? str ~ f3 n . Its contribution 
is small compared to that of j3 s , but it should not be overlooked for an unpaired 
system. For the oscillator potential, we have the following expression 

Jc 4 — 1 Oh 2 -I- 1 

+ J^ /3 +l AMR\ (16) 

It is seen that /3 n + /3 s tr < for the prolate nuclei with c/a < 3.15, whereas f3 8 is 
always positive. Thus, with an increase of the spin /, the ratio B/A has to change 
sign and to approach its limiting value B n /A n ~ A~ 8 ^ 3 , ~ 10 -6 for the SD bands 
in the 130 and 150 mass region, where rapid rotation destroys pairing correlations. 
The limiting ratio for a nucleus consisting of Z protons and iV neutrons is expressed 
by 



Bn (fc 4 -10fc 2 + l)P/ 3 

An ' (k 2 + 1) 3 A 8 / 3 



/2Zxl/3 /2iV\l/3 



(17) 
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One can therefore conclude that there are two distinct regions in the variation of 
B/A versus I. The lower part of a SD band is characterized by a gradual decrease 
of the pairing gap A. According to Eq. fll4|) the ratio B/A should exhibit a sharp 
increase. Then it changes sign and approaches the plateau ([T7|) at the top of a band 
because the deformation c/a depends weakly on spin in the normal phase. Such 
behavior of the B/A ratio is the signature of the pairing phase transition. 

We have analyzed all the SD bands of Ref. [|| with known or suggested spins 
of levels. Figure 2 shows the variation of the B/A ratio with / for bands with dif- 
ferent internal structure and different rotational frequencies. Apart from the bands 
192 Hg(l) and 194 Hg(3), where frequencies are so low that B/A rises continuously 
in the superffuid phase, and 84 Zr(l), for which pairing is quenched completely and 
B/A is close to the limiting value (|17D , all other bands display the behavior de- 
scribed above. It is important to note that such behavior is observed for the ND 
yrast band of 84 Zr, where B/A reaches the same limiting value (|17]) as in the SD 
band 84 Zr(l). There are other ND yrast bands of 168 Yb and 168 Hf with the phase 
transition which experimental evidence has been discussed previously in terms of the 
canonical variables |TT| and the spectrum of single-particle states |TT| , These 
bands exhibit the same features. Thus the plots of Fig. 2 demonstrate the uni- 
versality of the superffuid-to-normal phase transition for SD and ND bands. The 
manifestation of this universality was previously observed in the anomalous small 
value B/A — 7x 10~ 6 for the ND yrast band of 168 Hf fl3| and in the weak dependence 
of in the upper part of the SD band 152 Dy(l) |14j]. The new feature observed in 
this article is the small decrease of B/A at the top of the 152 Dy(l), 132 Ce(l), 84 Zr(l) 
and other SD bands. It may be explained by the decrease of the first multiplier in 
r.h.s. of Eq. (|17D due to increase of nuclear deformation at highest spins. 

All of the theoretical formulas obtained refer to collective rotation only. The 
influence of single-particle degrees of freedom on the B/A ratio is the subject of a 
separate investigation. Preliminary estimations show that the contribution of an odd 
nucleon in B s may be positive and irregular. Visible evidences of the nonrotational 
degree of freedom are shown in Fig. 3. A point of particular interest is the staggering 
band 149 Gd(l). It is seen that the staggering pattern is reproduced very well and the 
error bars are smaller than the staggering amplitude. We would like to emphasize 
that the physical values adequate for describing the staggering phenomenon are the 
B parameter or the B/A ratio. They do not require any smooth reference, which 
introduces an ambiguousness in the amplitude and the phase of oscillations. 

The next limit we want to consider is that of a large pairing gap A. In this 
case, the nonuniform pairing ( |T2"D is essential and the leading terms in the function 
$ are those proportional to the powers of Di and D\. They result in the limiting 
expression $ ~ (Huj /A) 2 . Thus for the very strong pairing (A ^> huj ), when the 
size of the Cooper pair Rhcu /A becomes much less than the nuclear radius, the 
value B s vanishes in agreement with the hydrodynamic equations of the ideal liquid 
|J. In the limit of an extremely large deformation, c/a — > oo, a needle shaped 
nucleus with pairing correlations rotates as a rigid body, S = Q r i g , B s = 0. For 
the finite but large deformation the deviations from these values are proportional to 
(a/c) 4 / 3 . This means that all nucleons with the exclusion of a small sphere in the 
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center of a nucleus are completely involved in rotational motion. Finally, for small 
deformations we have B s ~ (c/a — l) -6 that is comparable to the vibration-rotation 
interaction 0. 

Unlike the limiting value flTTD, it is impossible to compare with experiment the 
ratio B s /A s because the proton (A T ) and neutron (A„) pairing gaps are unknown 
for the SD bands. In such a case we try to solve an inverse problem. The equations 
for the two inertial parameters 

aA/% ig = Z<p(&, k) + N<pfa, k), (18) 

f3A/f3 = k) + k), (19) 

allow in principle to find £„. and £ w (^ = Ai/hu i, cu i = uj [2Ai/ A) 2 / 3 , I = ir,v). 
Here the function ip is taken from Ref. 



(p{v u u 2 ) = 1 - ? - r - — 2 2 — 2 , (20) 
(z/f + z^)[i/f^(i/i) + ^(^2)] 

and the rigid body moment of inertia and the value f3 Q have the form: 

% ig = 7.29A 5/3 {k 2 + l)k- 2/3 10- 3 h 2 MeV-\ 

P = 2.59A\k + l) 4 Ar 4/3 10- 8 ft 4 MeU~ 3 . (21) 
Finally the axis ratio k = c/a is determined by the quadrupole moment 

Qo = Q.05A 2/3 (k 2 - l)k- 2/3 l0~ 3 eb. (22) 



Unfortunately the system (18)— (19) does not have a solution. Table I gives the 
solutions of Eqs. flUf ) and ( |H| ) under simplifying assumption that the neutron and 
proton pairing gaps are equal. Different bands reveal the same tendency. The 
obtained values A x and A 2 show that the oscillator potential overestimates the 
moment of inertia and, to a greater extent, the (3 parameter. It is essential to use a 
more realistic potential for extracting pairing gaps from the inertial parameters. 



Table I. Pairing gap energies Ai and A2 (in MeV) found as the independent solutions 
of Eqs. (|T^ ) and (|l9| ) correspondingly in the case of ^ = The values c/a, ^s r i g , and f3 
were calculated by using Eqs. ([H|) and (|2^). The experimental values of a = and (3 were 
taken from Refs. § ( 172 Hf(yr)) and || ( 192 ' 4 Hg(l,2)). Those of 236 U(fi) were determined 
from Eqs. @, @. 



Band 


c/a %/% ig Ax p/fo A 2 


172 Hf(yr) 
192 Hg(l) 
194 Hg(l) 
236 U(fi) 


1.32 0.373 1.44 0.309 0.02 
1.61 0.762 0.95 0.151 0.05 
1.60 0.793 0.83 0.112 0.04 
1.84 0.816 1.01 0.141 0.08 
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In summary, the exact solution for the inertial parameter B in the superfTuid 
phase allows to show that neither super deformation nor nonuniform pairing arising 
from rotation induced pair density is responsible for the regularity of the SD ro- 
tational spectra. The extreme regularity of the SD bands in the 80, 130 and 150 
mass regions is explained by the transition from the superfluid to normal phase. 
The new signature of this transition reveals itself in the characteristic dependence 
of the ratio B/A with the spin I. Application of this criterion to experimental data 
indicates the existence of the phase transition in the SD bands of the three SD mass 
regions. A new extraction method of the proton and neutron pairing gap from the 
7-ray energies of interband transitions is discussed. 
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FIGURE CAPTIONS 



Fig. 1. Plot of the function $ from Eq. ([lj) against the dimensionless value £ for 
the axis ratio c/a = 2. The solid, dotted, and dashed lines correspond respectively 
to the exact value, the limit of close transitions, and the uniform pairing. The scale 
on the abscissa should be multiplied by a factor of approximately 7.7 for nuclei in 
the A ~ 150 mass region to obtain a gap energy in MeV. 

Fig. 2. B/A ratio versus spin for the SD and ND bands with mainly collective 
behavior. The formulas (|2|) are used to obtain this ratio from the experimental 
data for the ND (full circles) and SD (open circles) bands. The solid straight line 
is the limiting value B n /A n with the deformation c/a found from the quadrupole 
moment (p2|). Error bars (if they are greater than symbols) include 7-ray energy 
uncertainties only. Uncertainties in a spin assignment are not important for all the 
SD bands except 152 Dy(l), as the variation of spins in 2h would merely shift the 
curves along axis of abscissas. 

Fig. 3. Same as Fig. 2 for the bands with high-iV configurations having the 
non-zero alignment i. Straight line is the value of B n /A n for 153 Dy(l). 
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